Transition State Research Articles

Restriction on molecular fluxionality by substitution: A case study for the 1,10-dicyanobullvalene.

We show herein that 1,10-dicyano substitution restricts the paragon fluxionality of bullvalene to just 14 isomers which isomerize along a single cycle. The restricted fluxionality of 1,10-dicyanobullvalene (DCB) is investigated by means of: (i) Bonding analyses of the isomer structures using the adaptive natural density partitioning (AdNDP). (ii) Quantum dynamical simulations of the isomerizations along the cyclic intrinsic reaction coordinate of the potential energy surface (PES). The PES possesses 14 equivalent potential wells supporting 14 isomers which are separated by 14 equivalent potential barriers supporting 14 transition states. Accordingly, at low temperatures, DCB appears as a hindered molecular rotor, without any delocalization of the wavefunction in the 14 potential wells, without any nuclear spin isomers, and with completely negligible tunneling. These results are compared and found to differ from those for molecular boron rotors. (iii) Born-Oppenheimer molecular dynamics (BOMD) simulations of thermally activated isomerizations. (iv) Calculations of the rate constants in the frame of transition state theory (TST) with reasonable agreement achieved with the BOMD results. (v) Simulations of the equilibration dynamics using rate equations for the isomerizations with TST rate coefficients. Accordingly, in the long-time limit, isomerizations of the 14 isomers, each with Cs symmetry, approach the "14 Cs → C7v" thermally averaged structure. This is a superposition of the 14 equally populated isomer structures with an overall C7v symmetry. By extrapolation, the results for DCB yield working hypotheses for so far un-explored properties e.g. for the equilibration dynamics of C10H10.

Computational Analysis of Low Overpotential Ammonia Oxidation by Metal-Metal Bonded Ruthenium Catalysts, and Predictions for Related Osmium Catalysts.

The catalyzed electrochemical oxidation of ammonia to nitrogen (AOR) is an important fuel-cell half-reaction that underpins a future nitrogen-based energy economy. Our laboratory has reported spontaneous chemical and electrochemical oxidation of ammonia to dinitrogen via reaction of ammonia with the metal-metal bonded diruthenium complex Ru2(chp)4OTf (chp- = 2-chloro-6-hydroxypyridinate, TfO- = trifluoromethanesulfonate). This complex facilitates electrocatalytic ammonia oxidation at mild applied potentials of -255 mV vs ferrocene, which is the [Ru2(chp)4(NH3)]0/+ redox potential. We now report a comprehensive computational investigation of possible mechanisms for this reaction and electronic structure analysis of key intermediates therein. We extend this analysis to proposed second-generation electrocatalysts bearing structurally similar fhp and hmp (2-fluoro-6-hydroxypyridinate and 2-hydroxy-6-methylpyridinate, respectively) equatorial ligands, and we further expand this study from Ru2 to analogous Os2 cores. Predicted M24+/5+ redox potentials, which we expect to correlate with experimental AOR overpotential, depend strongly on the identity of the metal center, and to a lesser degree on the nature of the equatorial supporting ligand. Os2 complexes are easier to oxidize than analogous Ru2 complexes by ∼640 mV, on average. In contrast to mono-Ru catalysts, which oxidize ammonia via a rate-limiting activation of the strong N-H bond, we find lowest-energy reaction pathways for Ru2 and Os2 complexes that involve direct N-N bond formation onto electrophilic intermediates having terminal amido, imido, or nitrido groups. While transition state energies for Os2 complexes are high, those for Ru2 complexes are moderate and notably lower than those for mono-Ru complexes. We attribute these lower barriers to enhanced electrophilicity of the Ru2 intermediates, which is a consequence of their metal-metal bonded structure. Os2 intermediates are found to be, surprisingly, less electrophilic, and we suggest that Os2 complexes may require access to oxidation states higher than Os25+ in order to perform AOR at reasonable reaction rates.

Mechanism of Z-Selective Allylic Functionalization via Thianthrenium Salts.

A detailed mechanistic study of the Z-selective allylic functionalization via thianthrenium salts is presented. Kinetic analyses, deuterium labeling experiments, and computational methods are used to rationalize the observed reactivity and selectivity. We find that the reaction proceeds via a rate-determining and stereodetermining allylic deprotonation of an alkenylthianthrenium species. The Z-configuration of the resultant allylic ylide is translated into the Z-allylic amine product through a sequence of subsequent fast and irreversible steps: protonation to form a Z-allylic thianthrenium electrophile and then regioselective substitution by the nucleophile. In the stereodetermining deprotonation step, computational studies identified a series of stabilizing nonbonding interactions in the Z-alkene-forming transition state that contribute to the stereoselectivity.

Topological phase transition and tunable surface states in YBi.

A unique co-existence of extremely large magnetoresistance (XMR) and topological characteristics in non-magnetic rare-earth monopnictides has stimulated intensive research on these materials. Yttrium monobismuthide (YBi) has been reported to exhibit XMR up to 105% but its Topological properties still need clarification. Here we use the hybrid density functional theory to probe the structural, electronic, and topological properties of YBi in detail. We observe that YBi is topologically trivial semimetal at ambient pressure which is in accordance with reported experimental results. The topological phase transitions i.e., trivial to non-trivial are obtained with volumetric pressure of 6.5 GPa and 3% of epitaxial strain. This topological phase transitions are well within the structural phase transition of YBi (24.5 GPa). The topological non-trivial state is characterized by band inversions among Y-d band and Bi-p band at Γ- and X-point which is further verified with the help of surface band structure along (001) plane. The Z2 topological invariants are calculated with the help of product of parities and evolution of Wannier charge centers. The occurrence of non-trivial phase in YBi with a relatively small epitaxial strain, which a thin film geometry can naturally has, make it an ideal candidate to probe inter-relationship between XMR and non-trivial topology.

Spin-Restricted Descriptions of Singlet Oxygen Reactions from XMS-CASPT2 Benchmarks.

Reactions of singlet oxygen are numerous, some of which are desired but many are unwanted. Therefore, the ability to correctly predict and interpret this reactivity for complex molecular systems is essential to our understanding of singlet oxygen reactions. DFT is widely used for predicting many reactions but is not suited to degenerate electronic structures; application to isolated singlet oxygen often uses the spin-unrestricted formalism, which results in severe spin contamination. In this work, we demonstrate that spin-restricted DFT can correctly describe the reaction pathway for four prototypical singlet oxygen reactions. By careful benchmarking with XMS-CASPT2, we show that, from the first transition state onward, the degeneracy of the 1Δg state is broken due to differing interactions of the (degenerate) π* orbitals with the organic substrate; this result is well replicated with DFT. These findings demonstrate the utility of using spin-restricted DFT to explore reactions, opening the way to confidently use this computationally efficient method for molecular systems of medium to large organic molecules.

Compensación cinética y termodinámica del pardeamiento no enzimático de zumos clarificados de limón

This paper presents a study on kinetic and thermodynamic compensation to evaluate the kinetics of non-enzymatic browning of clarified lemon juices. In a previous work, the kinetic constants of browning were presented using the evolution of the absorbance at 420 nm (A420) and luminosity (L*) for different soluble solids contents (64.6; 50; 35; 20 and 10 ºBrix) and different working temperatures (70, 80, 90 and 95 ºC). The parameters of the Arrhenius equation were obtained by fitting the kinetic constants and temperatures. The variation of lnK0 with Ea follows a linear trend, so there is kinetic compensation, with isokinetic temperature values being 126.6 ºC and 149.7 ºC for A420 and L*, respectively. The transition state equilibrium constants were determined using the Eyring equation for each soluble solids content and each temperature. The equilibrium constants were adjusted to the Van't Hoff equation and the set of pairs of estimated values for the activation enthalpy and the activation entropy also followed a straight line, which causes thermodynamic compensation, with isoequilibrium temperatures of 112.8 ºC and 136.3 ºC for A420 and L*, respectively. It was concluded that the browning mechanism is the same for the ranges of soluble solids content and temperature studied. Since all isokinetic and isoequilibrium temperatures were higher than the working temperature values, it was also concluded that the control was enthalpic for all cases. To avoid the deterioration of these juices, it is advisable to intervene in the working temperature, trying to ensure that the heat treatment is carried out at the lowest possible temperature.

Combined Computational and Experimental Study Reveals Complex Mechanistic Landscape of Brønsted Acid-Catalyzed Silane-Dependent P═O Reduction.

The reaction mechanism of Brønsted acid-catalyzed silane-dependent P═O reduction has been elucidated through combined computational and experimental methods. Due to its remarkable chemo- and stereoselective nature, the Brønsted acid/silane reduction system has been widely employed in organophosphine-catalyzed transformations involving P(V)/P(III) redox cycle. However, the full mechanistic profile of this type of P═O reduction has yet to be clearly established to date. Supported by both DFT and experimental studies, our research reveals that the reaction likely proceeds through mechanisms other than the widely accepted "dual activation mode by silyl ester" or "acid-mediated direct P═O activation" mechanism. We propose that although the reduction mechanisms may vary with the substitution patterns of silane species, Brønsted acid generally activates the silane rather than the P═O group in transition structures. The proposed activation mode differs significantly from that associated with traditional Brønsted acid-catalyzed C═O reduction. The uniqueness of P═O reduction originates from the dominant Si/O═P orbital interactions in transition structures rather than the P/H-Si interactions. The comprehensive mechanistic landscape provided by us will serve as a guidance for the rational design and development of more efficient P═O reduction systems as well as novel organophosphine-catalyzed reactions involving P(V)/P(III) redox cycle.

Influence of Ti-incorporated Zeolite Topology and Pore Condensation on Vapor Phase Propylene Epoxidation Kinetics with Gaseous H2O2.

Vapor-phase propylene (C3H6) epoxidation kinetics with hydrogen peroxide (H2O2) strongly reflects the physical properties of Ti-incorporated zeolite catalysts and the presence of spectating molecules ("solvent") near active sites even without a bulk liquid phase. Steady-state turnover rates of C3H6 epoxidation and product selectivities vary by orders of magnitudes, depending on the zeolite silanol ((SiOH)x) density, pore topology (MFI, *BEA, FAU), and the quantity of condensed acetonitrile (CH3CN) molecules nearby active sites, under identical reaction mechanisms sharing activated H2O2 intermediates on Ti surfaces. Individual kinetic analyses for propylene oxide (PO) ring-opening, homogeneous diol oxidative cleavage, and homogeneous aldehyde oxidation reveal that secondary reaction kinetics following C3H6 epoxidation responds more sensitively to the changes in zeolite physical properties and pore condensation with CH3CN. Thus, higher PO selectivities achieved in hydrophilic Ti-MFI at steady-state reflect the preferential stabilization of transition states for C3H6 epoxidation (a primary reaction) relative to PO ring-opening and oxidative cleavage (secondary reactions) that solvation effects that reflect interactions among condensed CH3CN within pores and the extended pore structure.

Influence of Ti‐incorporated Zeolite Topology and Pore Condensation on Vapor Phase Propylene Epoxidation Kinetics with Gaseous H2O2

Vapor‐phase propylene (C3H6) epoxidation kinetics with hydrogen peroxide (H2O2) strongly reflects the physical properties of Ti‐incorporated zeolite catalysts and the presence of spectating molecules (“solvent”) near active sites even without a bulk liquid phase. Steady‐state turnover rates of C3H6 epoxidation and product selectivities vary by orders of magnitudes, depending on the zeolite silanol ((SiOH)x) density, pore topology (MFI, *BEA, FAU), and the quantity of condensed acetonitrile (CH3CN) molecules nearby active sites, under identical reaction mechanisms sharing activated H2O2 intermediates on Ti surfaces. Individual kinetic analyses for propylene oxide (PO) ring‐opening, homogeneous diol oxidative cleavage, and homogeneous aldehyde oxidation reveal that secondary reaction kinetics following C3H6 epoxidation responds more sensitively to the changes in zeolite physical properties and pore condensation with CH3CN. Thus, higher PO selectivities achieved in hydrophilic Ti‐MFI at steady‐state reflect the preferential stabilization of transition states for C3H6 epoxidation (a primary reaction) relative to PO ring‐opening and oxidative cleavage (secondary reactions) that solvation effects that reflect interactions among condensed CH3CN within pores and the extended pore structure.

Kinetic Computation of Cyclization Reactions of Large Keto-Hydroperoxide Radicals in Low Temperature Combustion.

The cyclization reactions of keto-hydroperoxide (KHP) radicals leading to the formation of keto cyclic ethers and OH radicals play an important role in low temperature combustion for hydrocarbon fuels or oxygenated hydrocarbon fuels. However, due to the lack of kinetic data of cyclization reactions of KHP radicals, researchers often derive high-pressure-limit rate constants of cyclization reactions of KHP radicals from analogous cyclization reactions of hydroperoxyl alkyl radicals during construction of the combustion mechanism. This study aims to systematically investigate the kinetics of cyclization reactions of KHP radicals involving short-to-large-sized radicals. The studied reactions are divided into 7 reaction classes, according to the size of the cyclic transition state, the conjugative effect (whether KHP radicals are resonance-stabilized or not), and the position of the carbonyl group (whether the carbonyl group is inside or outside of the reaction center). The isodesmic reaction method, in conjunction with transition state theory, is utilized for each reaction class to compute the energy barriers and high-pressure-limit rate constants at the DFT level. The study revealed that energy barriers calculated at the DFT level with correction by the isodesmic reaction method are close to the results from the benchmark CCSD(T) method. To develop more accurate rate rules, these reaction classes are further divided into subclasses based on the relative site of the OOH group with the carbonyl group, the type of carbon atoms where the OOH group is located, and the type of carbon atoms where the radical site is located. For each subclass, high-pressure-limit rate rules are derived by averaging the rate constants of reactions in the subclass, and it is found that the maximum absolute deviation of the energy barrier and the ratio of the largest rate constant to the smallest rate constant among reactions in each subclass are within chemical accuracy limits, indicating acceptable use of the developed rate rules. A comparison of the rate constants for cyclization reactions of KHP radicals with the values of analogous cyclization reactions of hydroperoxyalkyl radicals as provided in reported mechanisms is made. Additionally, a comparison is drawn between our developed rate rules for subclasses of the cyclization reactions of KHP radicals and the rate rules for analogous subclasses of cyclization reactions of hydroperoxyl alkyl radicals. These comparisons demonstrate significant differences and highlight the necessity for improved rate rules for cyclization reactions of KHP radicals to enhance the automatically generated combustion mechanisms for hydrocarbon and oxygenated hydrocarbon fuels.

Theoretical Studies on the Reaction Mechanism of Schiff Base Formation from Hexoses.

The Maillard reaction is one of the nonenzymatic post-translational modifications of proteins. Products of this reaction are considered to be related to aging diseases and the sensation of taste. In the initial stage of the Maillard reaction, Schiff base formation first occurs by the nucleophilic attack of amine nitrogen in proteins, and then, the reaction proceeds through the formation of 1,2-eminal and Amadori compounds. In this study, we computationally investigated the reaction pathway of Schiff base formation from hexoses. The optimized geometries of energy minima and transition states were calculated by using the density functional theory with the CAM-B3LYP/6-311+G(2d,2p) level of theory. The Schiff base formation progressed through three steps: two steps of carbinolamine formation and one step of dehydration. The dehydration is considered to be the rate-determining step in all hexoses because the activation barrier of the dehydration was higher than that of the carbinolamine formation. Furthermore, the steric configuration of the OH group at positions 2 and 3 affected the activation barrier.

Introduction of element Bi in the P-block promotes N2 activation for efficient electrocatalytic nitrogen reduction to produce ammonia

The synthesis of ammonia through electrocatalytic nitrogen reduction reaction (NRR) under ambient conditions is a presumable strategy. However, it still faces the problems of poor activity and low Faradic efficiency, so the development of efficient, stable and highly selective catalysts is of the importance. Molybdenum disulfide (MoS2) materials are widely used in NRR reactions owing to good electrical conductivity and abundant active sites, but the edge active site is the dual active site of NRR and competitive hydrogen evolution reaction (HER). Encouragingly, the strong interaction between elements in the P-block and N 2p orbitals results in excellent performance in adsorbing and activating nitrogen. Meanwhile, the partial occupation of p-orbitals by elements in the P-block also leads to poor HER activity. Herein, we introduced the P-block element Bi into 1 T-MoS2 and prepared Bi2S3/MoS2 composite. Taking advantage of the high conductivity of 1 T-MoS2 and the N2 adsorption and activation advantages of Bi element, an ammonia yield rate of 54.64 µg h−1 mg−1cat. and a high Faradic efficiency of 58.56 % were obtained. In addition, the active site and the evolution pathway of N2 reduction were investigated by in-situ infrared spectroscopy measurements and density functional theory (DFT). The in-situ infrared spectroscopy measurement results manifested the presence of the main intermediate transition state during the reduction of N2 to NH3. The calculation results indicated that N2 had the best adsorption on the introduced P-block element Bi, and the entire NRR reaction process followed an alternating pathway. This study reveals the synergistic effect of P-block elements and transition metal-based electrocatalytic materials on nitrogen reduction to ammonia synthesis.

Screening Stability, Thermochemistry, and Chemical Kinetics of 3-Hydroxybutanoic Acid as a Bifunctional Biodiesel Additive.

The thermo-kinetic aspects of 3-hydroxybutyric acid (3-HBA) pyrolysis in the gas phase were investigated using density functional theory (DFT), specifically the M06-2X theoretical level in conjunction with the cc-pVTZ basis set. The obtained data were compared with benchmark CBS-QB3 results. The degradation mechanism was divided into 16 pathways, comprising 6 complex fissions and 10 barrierless reactions. Energy profiles were calculated and supplemented with computations of rate coefficients and branching ratios over the temperature range of 600-1700 K at a pressure of 1 bar using transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) methods. Thermodynamics results indicated the presence of six stable conformers within a 4 kcal mol-1 energy range. The estimated chemical kinetics results suggested that TST and RRKM approaches are comparable, providing confidence in our calculations. The branching ratio analysis reveals that the dehydration reaction pathway leading to the formation of H2O and CH3CH═CHCO2H dominates entirely at T ≤ 650 K. At these temperatures, there is a minor contribution from the simple homolytic bond fission reaction, yielding related radicals [CH3•CHOH + •CH2CO2H]. However, at T ≥ 700 K, this reaction becomes the primary decomposition route. At T = 1700 K, there is a minor involvement of a reaction pathway resulting in the formation of CH3CH(OH)•CH2 + •CHO(OH) with an approximate contribution of 16%, and a reaction leading to [•CH3 + •CH2OHCH2CO2H] with around 9%.

Darobactin Substrate Engineering and Computation Show Radical Stability Governs Ether versus C-C Bond Formation.

The Gram-negative selective antibiotic darobactin A has attracted interest owing to its intriguing fused bicyclic structure and unique targeting of the outer membrane protein BamA. Darobactin, a ribosomally synthesized and post-translationally modified peptide (RiPP), is produced by a radical S-adenosyl methionine (rSAM)-dependent enzyme (DarE) and contains one ether and one C-C cross-link. Herein, we analyze the substrate tolerance of DarE and describe an underlying catalytic principle of the enzyme. These efforts produced 51 enzymatically modified darobactin variants, revealing that DarE can install the ether and C-C cross-links independently and in different locations on the substrate. Notable variants with fused bicyclic structures were characterized, including darobactin W3Y, with a non-Trp residue at the twice-modified central position, and darobactin K5F, which displays a fused diether ring pattern. While lacking antibiotic activity, quantum mechanical modeling of darobactins W3Y and K5F aided in the elucidation of the requisite features for high-affinity BamA engagement. We also provide experimental evidence for β-oxo modification, which adds support for a proposed DarE mechanism. Based on these results, ether and C-C cross-link formation was investigated computationally, and it was determined that more stable and longer-lived aromatic Cβ radicals correlated with ether formation. Further, molecular docking and transition state structures based on high-level quantum mechanical calculations support the different indole connectivity observed for ether (Trp-C7) and C-C (Trp-C6) cross-links. Finally, mutational analysis and protein structural predictions identified substrate residues that govern engagement to DarE. Our work informs on darobactin scaffold engineering and further unveils the underlying principles of rSAM catalysis.

Photoionization and Dissociative Photoionization of Acetaldehyde in the 10.0-13.7 eV Range by Synchrotron Photoelectron Photoion Coincidence Spectroscopy.

Photoionization and dissociative photoionization of acetaldehyde (CH3CHO) in the 10.0-13.7 eV energy range are studied by using synchrotron radiation double imaging photoelectron photoion coincidence spectroscopy (i2PEPICO). The X2A' and A2A" electronic states of CH3CHO+ as well as the Franck-Condon gap region between these two states have been populated with several vibrational sequences and assigned in the high-resolution slow photoelectron spectrum (SPES). The adiabatic ionization energies (AIEs) of the X2A' and A2A" states are measured at 10.228±0.006 and 12.52±0.05 eV, respectively. The present results show that the X2A' state is a stable state while the A2A" state is fully dissociative to produce CH3CO+, CHO+ and CH4 + fragment ions. The 0 K appearance energies (AE0K) of CH3CO+ and CHO+ fragment ions are determined through the modeling of the breakdown diagram, i. e., AE0K(CH3CO+)=10.89±0.01 eV (including a reverse barrier of ~0.19 eV) and AE0K(CHO+)=11.54±0.05 eV. In addition, the dissociation mechanisms of CH3CHO+ including statistical dissociation, direct bond breaking and isomerization are discussed with the support of the calculated dissociation limits and transition state energies.

Effects of sulfur‐ and Fe‐deficient defects on the decomposition of hydrogen peroxide H2O2 on pyrite surface

AbstractThere is still a lack of research on the effect of the S and Fe vacancy on the adsorption mechanism of peroxide (H2O2) on the pyrite surface. In this work, the adsorption of H2O2 molecules on the perfect and defective pyrite surface is simulated based on density functional theory (DFT). The results indicated that the adsorption energy of H2O2 molecules on the S‐vacancy surface is relatively lower than that on the perfect surface and dissociative adsorption entirely occurs, while the adsorption behavior of H2O2 molecules on the S‐rich pyrite surface has not occurred. When the H2O2 molecules are adsorbed on the Fe‐rich pyrite surface, two products appear, one with two OH− species and the other with an H2O molecule and an O atom. The analysis of the transition state (TS) found that the dissociation of the H2O2 molecule producing two OH− species that bond with two Fe atoms is more likely to occur. In addition, based on the coordination chemistry, it is demonstrated that Fe atoms are more reactive around an S vacancy than on a perfect surface.

Experiment and Computational Study on Pd‐Catalyzed MIA‐Directed Ortho‐C‐H Alkynlation of Phenylalanine

The direct ortho‐alkynylation of substituted phenylalanine has been successfully achieved through methoxyiminoacyl (MIA)‐mediated Pd‐catalyzed C‐H functionalization. This innovative protocol showcases the ability to apply a diverse range of phenylalanine substrates, resulting in the efficient synthesis of alkynylation benzylamine derivatives with notable effectiveness. Computational investigations further revealed that the fluoride ion enhances the electropositivity of Pd(IV) and the interaction between C1 and C2 within the transition state of the reductive elimination, which significantly expedites the reductive elimination step in the alkynylation process.

Turning noncovalent interactions via halogen substituents for improving activity and enantioselectivity of Candida antarctica lipase B

Noncovalent interactions associated with the enzyme, enzyme-substrate complexes, and transition states have profound influences on varying the biocatalytical activity and enantioselectivity. The recent rediscovery of halogen bonds incorporating in biocatalysis has induced our motivation on introducing halogen substituents into the substrate for improving the enzyme performance. By selecting Candida antarctica lipase B (CALB)-catalyzed kinetic hydrolysis (or butanolysis) of trans-2-(p-halophenyl)cyclopropyl 1,2,4-triazolide in methyl tert-butyl ether (MTBE) as the model reaction system, (R, R)-2-(p-halophenyl)cyclopropane-1-carboxylic acid (or ester) of high-optical purity that acts as the key building block for synthesizing the corresponding (1 R, 2S)-2-(p-halophenyl)cyclopropylamine derivatives was obtainable. The thermodynamic and kinetic analysis was then employed for elucidating effects of the halogen substituent and acyl acceptor on enhancing the noncovalent interactions of the rate-limiting tetrahedral intermediate and hence CALB activity and enantioselectivity. A putative model was moreover proposed for interpreting the possible hydrogen and halogen bonding involving the halogen substituent in turning the noncovalent interactions.

Probing the activated complex of the F + NH3 reaction via a dipole-bound state

Experimental characterization of the transition state poses a significant challenge due to its fleeting nature. Negative ion photodetachment offers a unique tool for probing transition states and their vicinity. However, this approach is usually limited to Franck-Condon regions. For example, high-lying Feshbach resonances with an excited HF stretching mode (vHF = 2-4) were recently identified in the transition-state region of the F + NH3 → HF + NH2 reaction through photo-detaching FNH3− anions, but the direct photodetachment failed to observe the lower-lying vHF = 0,1 resonances and bound states due apparently to negligible Franck-Condon factors. Indeed, these weak transitions can be resonantly enhanced via a dipole-bound state (DBS) formed between an electron and the polar FNH3 species. In this study, we unveil a series of Feshbach resonances and bound states along the F + NH3 reaction path via a DBS by combining high-resolution photoelectron spectroscopy with high-level quantum dynamical computations. This study presents an approach for probing the activated complex in a reaction by negative ion photodetachment through a DBS.

Physical Prior Mean Function-Driven Gaussian Processes Search for Minimum-Energy Reaction Paths with a Climbing-Image Nudged Elastic Band: A General Method for Gas-Phase, Interfacial, and Bulk-Phase Reactions.

The climbing-image nudged elastic band (CI-NEB) method serves as an indispensable tool for computational chemists, offering insight into minimum-energy reaction paths (MEPs) by delineating both transition states (TSs) and intermediate nonstationary structures along reaction coordinates. However, executing CI-NEB calculations for reactions with extensive reaction coordinate spans necessitates a large number of images to ensure a reliable convergence of the MEPs and TS structures, presenting a computationally demanding optimization challenge, even with mildly costly electronic-structure methods. In this study, we advocate for the utilization of physically inspired prior mean function-based Gaussian processes (GPs) to expedite MEP exploration and TS optimization via the CI-NEB method. By incorporating reliable prior physical approximations into potential energy surface (PES) modeling, we demonstrate enhanced efficiency in multidimensional CI-NEB optimization with surrogate-based optimizers. Our physically informed GP approach not only outperforms traditional nonsurrogate-based optimizers in optimization efficiency but also on-the-fly learns the reaction path valley during optimization, culminating in significant advancements. The surrogate PES derived from our optimization exhibits high accuracy compared to true PES references, aligning with our emphasis on leveraging reliable physical priors for robust and efficient posterior mean learning in GPs. Through a systematic benchmark study encompassing various reaction pathways, including gas-phase, bulk-phase, and interfacial/surface reactions, our physical GPs consistently demonstrate superior efficiency and reliability. For instance, they outperform the popular fast inertial relaxation engine optimizer by approximately a factor of 10, showcasing their versatility and efficacy in exploring reaction mechanisms and surface reaction PESs.